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2^ ■ Abstract 

T-H ' The parametrization of the reaUstic Bethe-Salpeter ampUtude for the deuteron 

is given. Eight components of the ampHtude in the Eudidean space are presented 
^ . as an analytical fit to the numerical solution of the Bethe-Salpeter equation in the 

^ I ladder approximation. An applicability of the parametrization to the observables of 

^ ■ the deuteron is briefly discussed. 

I ! 

^ ■ 1 Motivations 

(N : 

0> I Since the first solution of the realistic Bethe-Salpeter (BS) equation for the deuteron the 

ly-^ ■ BS-amplitudes have been applied to describe various processes with the deuterons |]I|, 0, H, 

O ■ Hi- The obvious advantages of the approaches based on the BS formalism are the explicit 

^ ! covariance and connection to the covariant dynamical (field) theory. In spite of this, practical 

use of the BS amplitudes is not as popular as of nonrelativistic wave functions [P, A more 
Q^l complicated physical interpretation and technical complexity of the approaches based on the 

p • BS amplitudes are the main reasons for that. 

^ ■ In a series of recent papers it has been argued that a new intuition can be developed 

• • . in working with the BS amplitudes 0, |^. It has been shown that calculations of many 

i observables for the deuteron are reduced to calculations similar to those in field theory. It 

^ I has also been stressed that the usage of the BS amplitude for the deuteron in the Dirac 

' matrix basis can be more convenient than in the spinor basis. In this case computations are 

formalized enough to extensively apply analytical computing software, such as the Mathe- 
matica [Q, to calculate the matrix elements of observables in terms of the components of 
the BS amplitude. However, the promotion of a new technique should assume that all ba- 
sic ingredients of the calculation are available to a potential user. The BS formalism still 
lacks this feature, since there is no simple parametrization available for the realistic deuteron 
amplitude, such as is available for nonrelativistic wave functions [0, ^ or relativistic wave 
functions of the spectator equation |[TI|. 



This letter presents the analytical parametrization of the Bethe-Salpeter amplitude for 
the deuteron. The parameters are fixed by fitting to the recent numerical solution of the 
homogeneous BS equation using the Dirac matrix basis 0. The one-boson exchange 
potential from ref. |l^ was used with a minor adjustment of its parameters 0, so in this 
sense the solution does not contain new physics (or different physics) then the pioneer paper. 
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2 Definitions and kinematics 



The realistic BS amplitude of the deuteron, $, can be obtained as a solution to the homo- 
geneous BS equation with the effective one-boson-exchange kernel |jl|, |^: 

4(p. P.) ^ .Ste.p.) E / (^ (/>_^^.^., ^(P'. i-.). (1) 

where Po is the deuteron momentum, /is is the mass of the meson 5, F^'' is the meson- 
nucleon vertex, corresponding to the meson B and connected to the fc-th nucleon. The tensor 
product in the r.h.s. of eq. (|^) is for a possible interconnection of the quantum numbers 
in two vertices (such as in the case of the vector or isovector mesons). The two-nucleon 
propagator, 5(^1,^2), is defined as: 

•5(^1,^2) = 2 ^^^2"! • = {p^ + m){p2 + m)D{p, Pd), (2) 

Pi — + le P2 — ni + le 

Pl,2 = ^±P, (3) 

where m is the nucleon mass, pi^2 are nucleon momenta, pk = l^iPk -^(P) Pd) is "scalar 
two particle propagator". 

It is convenient to use the amplitude, \E', conjugated with respect to one of the nucleon 
lines §: 

$ = vl>7^ or \E' = -$7^ (4) 

where 7^^ = 7371. Then the amplitude \E' can be decomposed in terms of a complete set of 
the Dirac matrices, their bilinear combinations and the 4 x 4-identity matrix, 1, with sixteen 
components being the coefficients of decomposition: 

^ = ^2 + l^^v + + l^li^Va + C^/..V^r- (5) 

In the deuteron rest frame the notions are used: 

C^lCV'J, C = (€,V'J, (6) 

where j. A; = 1, 2, 3 and other tensor components of V'f''' are equal to zero. 

In order to separate the amplitude with the deuteron's quantum numbers, a partial wave 
decomposition of the four vector and four scalar functions, (|^)-®, is performed: 

^= EV-IpcIpI; J)Yjm(^^p) (8) 

JM 
JM 

+ v^(po, IpI; JJ)Y'j^,{<:i,) + 7A(po, IpI; J + \J)y'jX}{<:i,)\ , (9) 



where Yjm and Yjm are the spherical harmonics and vector spherical harmonics respectively. 

Fixing the total momentum J = 1 and separating the components with positive parity, 
the deuteron's amplitude components read as: 

^O(l), ^,(11), ^,(01), ^,(21), tp^ill), ^,(01), M^l). (10) 

The components ipa{l) and ip^lll) are odd functions of po and all others are even. 

The BS equation with the realistic one-boson exchange potential is solved using the Wick 
rotation, which corresponds to replace po — *• ip^ and V'a(l) ~^ ^V'a(l) O- This procedure 
removes singularities from the exchange meson propagators in eq. (|lD and from the scalar 
propagator, D, which in the deuteron rest frame takes the form: 

Dip,Pn)-' = Dip,, \p\)-' = [m' + p'+pl - ^Mlf + plMl, (11) 

where Mr, = 2m + e^) is the deuteron mass. 

After the Wick rotation, the components of the deuteron amplitude are computed along 
the imaginary axe in the complex po-plane. Since the inverse Wick rotation of the numeri- 
cally known amplitude is an ill-defined operation, the parametrization for the components is 
obtained in the Wick rotated case. The possibility to analytically continue those amplitudes 
into a physical region is discussed in Section 4. 



3 The parametrization 

The parametrization of all components has the form (index J = 1 is omitted): 

^(P4, \p\]L) = f{p4, \p\;L)Exp [g{\p\; L)pI] D{p,, |^'|)m^ (12) 

where L = J, J ±1, depending on the quantum numbers of the component and functions / 
and g are given by: 

I IL , ,L 

f{p,Ap\;L)=j:MP^)-^,, 9mL) = Y.B^-^fl—,. (13) 

The form of parametrization (p!3|) and the scale of the parameters ai are prompted by the 
previous works with parametrizations of the wave functions |^, [1^, |Tl|: 

ao = /io/V2, ai = i^o, i = l...Nf, 

Pi = ifii, i = l...Ng, 

Ho = 0.139 GeV, /ii = 2/io, 

m = 0.939 GeV, Md = 2m + en, to = -2.2246 MeV. (14) 

The number Nf is equal to 11 for all components, whereas Ng differs for different components. 
The p4-dependence of the coefficients Ai^p,) is given by: 

AM=A+pIA',, for V'p(l),^.(l),^a(0),^a(2),^t(0),^i(2); 

A{p,)=p,iA+plA'^, for V^°(1),V^?(1). (15) 
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Coefficients Ai, A[ and Bi for all components are given in Appendix A, Tables |T1-|. The 
presented parametrization contains a seemingly large number of parameters, 27 to 29 per 
every of eight components plus three parameters common for all of them, including the 
nucleon mass, m, deuteron binding energy, e^j, and an additional mass scale parameter, 
/io- This looks rather unusual for such type of parametrizations. The parametrizations of 
wave functions |]^, [1^, |TT|, for instance, contain only n ~ 10 parameters per component. 
The reason for this difference is that the components of the BS amplitude depend upon two 
independent variables, the relative momentum, p, and "relative energy", p^. It is clear now, 
that the presented parametrization contain quite a modest number of parameters; it is not 
even compared to the one dimensional fit of the wave functions. 



4 The applicability of the parametrization 



The parametrization (|12D-(|15]) is obtained by fitting the numerical solution to the BS equa- 
tion, using the least-squares procedure The domain of validity of the parametrization 
in relative momentum is < |p| < 3 GeV, means that the solution of the BS equation was 
fitted up to this point. The domain of validity in relative energy (which is actually ipo) is 
defined as follows. First, the singular structure of the BS amplitude in the Minkowski space 
is governed by the singularities of the propagator D{Pr),p), eqs. (g) and (|11|), where the 



closest nucleon pole is most important for the physical applications. Thus, the parametriza- 



tion is valid at least up to p4 ~ Md/2 — ^m? + corresponding to the nucleon pole at 
given IpI . Second, the parametrization allows for the integration in the matrix elements over 
Pi with infinite limits, (—oo, +cxd), in the Euclidean space. 

The starting point for calculating any quantity with the BS amplitude is the relativistic 
impulse approximation. In many cases the relativistic impulse approximation is presented 
by the Feynman "triangle diagram" with zero transfer of the momentum ||^, ^ ^, |^, g = 
(Fig. P: 

(O) = / ^Tr{^'(po,p)0vl>(po,p)(p2-m)}, (16) 

where ^ = 7o\&^7o. Two important examples are the matrix elements of the vector and axial 
currents, 6 = 7^, 757^,. 

The matrix element (70), the vector charge, is used to normalize the BS amplitude: 



2Mn I (27r)4 

dp4d\p\p'^ 



1 



Tr{^'(po,p)7o*(po,p)(p2 -m)} (17) 

{ 



-8m {Mo)Mo)+M'^)M'^)) 
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+(M, - 2p,) (v^°(i)' + Mof + M^f + M^f 

+Mlf + AMOf + 4^t(2)' + 4^?(1)') } ■ (18) 

The components in eq. (plsj) are parametrized by eqs. (p!^)- ([l5|) . (Note that factor 2tt from 
integration over angle is absorbed by the parametrization.) Integrating over in eq. (|18|) 
but keeping |p|-dependence, one gets the charge density in the momentum space, analogous 
to the square of the deuteron wave function in a nonrelativistic approach. This charge 
density is shown in Fig. ^ together with the density of the 3-rd component of the axial 
current (omitting terms, vanishing after integration over 6): 

(7573)= 2]fe /^^Tr{^'(po,p)7573^(po,p)(P2 -m)} (19) 



dp4d\p\p 



1 



2 

-{ 



(2vr; 

2m (4^,(0)^t(0) - 2^,(2)^,(2) 

-2v^V^°(l)7/;°(l) + V2^,{1)M^: 

+ ^ (4^p(l)^,(0) + 2v^^p(l)^,(2) 



-2v^^,(0)^,(l)+ ^,(2)^,(1) 

+ (M, - 2p4) (-V^a(O)' + ^Va(2)' - ^M^? 

-AMOf + m^f-^i^^af)]- (20) 

The quality of the parametrization has been checked by a comparison of the charge and 
axial densities, as well as their integrals, computed using the parametrization and original 



numerical components. The original amplitude is normalized by (^Tj) and is exactly equal 
to 1.0, whereas the parametrization gives the normalization equal to 0.9997. This is not a 
trivial result, since all components were fitted independently. One can use this number for 
the "renormalization" of observables. The original amplitude gives an axial charge value 
of 0.9215, while the parametrization yields the same. The error of the parametrization 
describing the densities is ~ 0.01% at small |p| to ~ 1 — 2% at |p| ~ 1-3 GeV. 

Finally, the issue of an "inverse Wick rotation" should be addressed. The analysis of 
the singular structure of the "triangle graph" and behavior of the BS amplitude result in 
the conclusion that, perhaps, the presented parametrization can be used for an analytical 
continuation, p^ —ipo, of the BS amplitude up to the closest nucleon pole of po = Md/2 — 



m? + However, such a procedure works well (with accuracy ~ 10%) only up to 

IpI ~ m. The accuracy was estimated by calculating the vector and axial densities for the 
processes with the second nucleon on mass-shell. If it is used further, the procedure gives 
an accuracy of ~ 50% at |p| = 1.5 GeV and it should not be used beyond this point. 



5 



5 Summary 

The parametrization of the reahstic Bethe-Salpeter amphtude for the deuteron has been 
presented. All eight components of the amplitude are given in the Wick rotated case in the 
form of analytical functions. Simple examples of the use of the parametrization are presented 
and the applicability domain is discussed. 
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Appendix A. Tables of parameters 

Table 1: Parameters for the V'p(l) component 



L 


i 


A 




Bi 


1 





-0.0108242476 


-34.3742135 






1 


0.0576888290 


119.943051 







2 


-3.00806368 


-1257.32744 







3 


52.5918636 


13865.0551 


-222.089494 




4 


-549.793133 


-101125.416 


890.102891 




5 


3325.75711 


472699.731 


-1210.29749 




6 


-11944.6774 


-1417223.72 


548.300818 




7 


26383.0580 


2743193.30 






8 


-34689.7448 


-3406750.03 






9 


24977.5913 


2621309.83 






10 


-8159.59409 


-1137897.83 






11 


607.051002 


213102.794 





Table 2: Parameters for the ^a(l) component 



L 


i 


A, 


A[ 


B^ 


1 





0.113315679 


15.2791178 






1 


-0.565717405 


-42.6551270 







2 


12.3795242 


292.530922 







3 


-211.521767 


-2058.70255 







4 


2025.71260 


8992.57139 


-48.3735989 




5 


-11769.4095 


-20393.3149 


91.9530671 




6 


43037.5609 


10486.7684 


-44.7438523 




7 


-101273.551 


55592.5995 






8 


151268.732 


-141733.721 






9 


-136889.618 


152338.654 






10 


67943.7391 


-80778.8808 






11 


-14143.4907 


17289.2804 
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Table 3: Parameters for the ipvi^) component 



L 


i 




A', 

t 


Bi 


1 





00545480918 


452976149 






1 


-0 0413469671 


3 58498830 







2 


-0 626170875 


49 3118036 







3 


-12 5904324 


-435 917083 







4 


144.575745 


3114.24317 


6.14953994 




5 


-1128.73423 


-14069.8986 


-25.2161192 




6 


4901.27381 


41906.4275 


19.9844772 




7 


-13706.4378 


-85174.0821 






8 


24829.8901 


114516.487 






9 


-26877.9006 


-95556.5511 






10 


15516.3123 


44352.4015 






11 


-3665.63223 


-8706.31030 





Table 4: Parameters for the V^a(O) component 



L 


i 


A^ 


A[ 


5^ 








-0.00198690442 


-2.91350897 






1 


0.01915782253 


12.2001594 


-0.853805793 




2 


-0.154455300 


-214.362824 


17.7585282 




3 


25.2282913 


3289.28876 


-89.3945864 




4 


-361.376166 


-30582.8134 


155.578411 




5 


2946.98165 


166962.671 


-92.0411385 




6 


-13477.6590 


-553839.568 






7 


37613.2426 


1152552.96 






8 


-63648.9139 


-1515457.41 






9 


62158.0044 


1223772.96 






10 


-31778.3109 


-554471.928 






11 


6523.12231 


107982.421 
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Table 5: Parameters for the V'a(2) component 



L 


i 




A', 


Bi 


2 





0906145663 


-54 7995054 






1 


-0 712220945 


256 049408 


-57 8195768 




2 


-20 3136219 


470 222091 


211 630253 




3 


81 4188182 


-8294 9022 


-385 318964 




4 


-544.595070 


33870.9853 


367.461861 




5 


2355.01212 


-77139.6959 


-137.187405 




6 


-6864.79033 


118585.336 






7 


13318.8325 


-138016.983 1 






8 


-16303.3461 


126552.418 






9 


11929.7419 


-84889.6419 






10 


-4708.14062 


35162.7881 






11 


756.740280 


-6501.81818 





Table 6: Parameters for the V'°(l) component 



L 


i 


Ai 


A'i 


Bi 


1 





0.0126177910 


0.547674806 






1 


-0.0622969395 


-3.75788455 







2 


1.64026310 


30.8407337 







3 


-25.8788379 


-220.882836 







4 


241.662795 


1370.42460 


7.63938144 




5 


-1410.95190 


-5648.39611 


-31.1836535 




6 


5126.52517 


15857.5305 


24.3738723 




7 


-12077.1617 


-30611.8247 






8 


18231.6691 


39147.3503 






9 


-16701.2532 


-31118.3495 






10 


8364.33221 


13774.6556 






11 


-1750.41655 


-2578.09116 





9 



Table 7: Parameters for the ■0t(O) component 



L 


i 


A., 


A'i 


Bi 








-0 000764090007 


-0 909552636 






1 


00774111957 


3 53114793 


0266757851 




2 


081210854 


-79 1547771 


4 16947945 




3 


8 49620518 


1058 66511 


-3 26685212 




4 


-122.622299 


-8772.66567 


-36.1206906 




5 


1056.75816 


43565.4328 


36.1389858 




6 


-5058.22423 


-132777.961 






7 


15091.7728 


258802.623 






8 


-28171.6682 


-327964.927 






9 


31347.0373 


263416.014 






10 


-18926.2139 


-122064.787 






11 


4783.25176 


24814.8930 





Table 8: Parameters for the V^t(2) component 



L 


i 


Ai 




5^ 


2 





0.0927963894 


6.63409360 






1 


-0.566806578 


27.0403716 


-26.4480841 




2 


-6.23310249 


628.891454 


39.0979395 




3 


-22.0317803 


-5172.06340 


-12.4959784 




4 


299.860807 


22767.0450 






5 


-2053.92122 


-74280.5598 






6 


8124.81954 


177690.047 






7 


-20191.0789 


-294391.655 






8 


31928.1425 


323661.633 






9 


-30700.6358 


-224618.740 






10 


16220.4293 


89089.3349 






11 


-3598.88935 


-15407.7777 
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Figure captions 

Figure 1: The diagram for the matrix element of operator O over the deuteron state in the 
impulse approximation. 



Figure 2: Densities of the vector (sohd hne) and axial (dashed line) charges calculated with 
the presented parametrization. 
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Fig. |I[ A. Umnikov. 
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Fig. ^ A. Umnikov. 
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